A dynamical collective calculation of supernova neutrino signals 
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We present the first calculations with three flavors of collective and shock wave effects for neutrino 
propagation in core-collapse supernovae using hydroynamical density profiles and the S matrix 
formalism. We explore the interplay between the neutrino-neutrino interaction and the effects of 
multiple resonances upon the time signal of positrons in supernova observatories. A specific signature 
is found for the inverted hierarchy and a large third neutrino mixing angle and we predict, in this 
case, a dearth of lower energy positrons in Cherenkov detectors midway through the neutrino signal 
and the simultaneous revelation of valuable information about the original fluxes. We show that this 
feature is also observable with current generation neutrino detectors at the level of several sigmas. 
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Introduction Core-collapse supernovae are rare but rep- 
resent impressive violent phenomena in the Universe ca- 
pable of emitting 10^^ ergs as neutrinos of all flavors in 
a short burst of about 10 seconds. While such neutri- 
nos were first observed in 1987, current observatories or 
next generation detectors, could provide us with the lu- 
minosity curve from a future (extra)galactic explosion 
and/or the exciting observation of relic neutrinos from 
past supernovae. Untangling the information from such 
neutrino signals will represent a challenging task, since 
the information from the explosion phenomenon and the 
neutrino properties - e.g. the neutrino hierarchy and the 
third neutrino mixing angle - is intertwined. 

Impressive progress has been achieved in the last few 
years in our understanding of how neutrinos propagate in 
supernovae, fundamentally modifying the standard MSW 
effect paradigm. This evolution is due to the substantial 
progress made in the calculations which now either in- 
clude neutrino-neutrino interactions or use profiles which 
contain multiple resonances. Pantaleone first pointed out 
that the neutrino-neutrino interactions lead to an off- 
diagonal coupling [Ij] , and subsequent numerical calcula- 
tions showed that new phenomena arose when such terms 
are included in the neutrino evolution equations 01- It 
is now well established that neutrinos evolve collectively 
through the synchronized and bipolar regimes eventu- 
ally exhibiting a spectral split (see e.g. 0, 0, H, @, 0]). 
The phcnomcnological implications are important, espe- 
cially for the inverted hierarchy case, because electron 
anti-neutrinos get a "hot" spectrum, independently of the 
smallness of ^13 . A second swapping of the anti-neutrino 
fluxes may also occur when they reach the MSW reso- 
nance region, depending on the hierarchy and of the size 
of . Recent works have pointed out that the dynamics 
of the explosion can alter the high density resonance Q 
and have found that the neutrinos might encounter the 



resonance density more than once (multiple resonances) 
13, 111- Rapid oscillations in the (anti)neutrino propa- 
bilities might appear because of the phase differences be- 
tween two or more resonances [ill, [l3| • 

This paper presents the first numerical calculation with 
three flavors to include both the neutrino-neutrino in- 
teraction and dynamic MSW effects using matched den- 
sity proflles and correctly putting them together using 
S matrices rather than probabilities. We focus upon 
the results of the anti-neutrino time signal. Signatures 
pinpointing to the hierarchy and possible ^13 values are 
shown, even for values beyond the proposed reach of fu- 
ture experiments on Earth. For large ^13 and inverted 
hierarchy we predict that the positron event rate pos- 
sesses a characteristic feature due to the passage of the 
shock. Such a feature is observable in current and future 
large-size neutrino observatories. 

Theoretical framework Wc calculate the three flavor neu- 
trino evolution through matter in two steps. First, we 
evolve the neutrino wavefunctions up to some radius us- 
ing supernova density profiles at different times during 
the supernova explosion, as done in [lit for one static 
density profile. This calculation includes the neutrino 
coupling both to matter with loop corrections (i.e. V^t) 
and to neutrino themselves. For the latter we use the 
single-angle approximation, i.e. wc assume that neu- 
trinos are essentially emitted with one angle. Such an 
assumption accounts rather well both qualitatively and 



quantitatively for the neutrino collective effects la, 
even though in some cases decoherence in a full multi- 
angle description might appear (see e.g. [13 )• The den- 
sity profile used is a dynamic inverse power-law. The 
second step is to determine the exact neutrino evolution 
through the rest of the supernova mantle by solving the 
evolution operator equations as described in [l3| which 
is a 3 flavor generalization of 11[. The ID density pro- 
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FIG. 1; Antineutrino oscillation probabilities, as a function 
of the distance from the neutron-star surface, including the 
vv interaction and the V^t- The curves correspond to anti- 
neutrinos having 15 MeV energy, the hierarchy is inverted and 
^13 is large. 



files used are taken from the hydrodynamical calculations 
described in 1J| and include both a front and a reverse 
shock. These profiles are matched to the profiles used in 
the first step. The full results are then sewn together us- 
ing the S matrices rather than probabilities c.f. 
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Finally, the flux on Earth is calculated taking into ac- 
count decoherence [l^l but not Earth matter 21 1 which 
would be relevant if the supernova were shadowed. 

Numerical results Our main goal is to explore the neu- 
trino time signal in an observatory, depending on the 
yet unknown neutrino parameters, and see if we can ex- 
ploit a combination of the neutrino-neutrino interaction 
and shock wave effects to get clues on important open 
issues. We take as an example the electron anti-neutrino 
scattering on protons which is the dominant channel in 
Chercnkov and scintillator detectors. The results we 
present are obtained with the best fit oscillation param- 
eters, i.e. Ato?2 = 8 X lO-^eV^, sin^ 26*12 = 0.83 and 
jAmlgl 3 X 10~^cV2, sin226'23 = 1 for the solar and 
atmospheric differences of the mass squares and mixings, 
respectively [2^. The Dirac phase is taken to be zero: 
no effects show up when the muon and tau luminosities 
are taken equal |23| ; while a few percent modification can 
appear due to and of non- linear effects [15 1. 

One of the important open questions is the hierarchy, 
since the sign of Arri^^ can be positive (normal) or nega- 
tive (inverted hierarchy) . The value of the third neutrino 
mixing angle is another issue, particularly crucial for the 
search of CP violation in the lepton sector. We take here 
two possible values for ^13, a large (sin^^ia = 10~^) and 
a small (sin^^is = 10~®). Note that the results corre- 
sponding to the large value are emblematic of the whole 
range sin'^^ia > 10~^ up to the present experimental 
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FIG. 2: Electron anti-neutrino fluxes on Earth, in the case of 
adiabatic ("cold" spectrum, solid) and non-adiabatic ("hot" 
spectrum, dashed) conversions in the star. The time is 1 s. 
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FIG. 3: Positron time signal per unit tonne in a detector, for 
a galactic explosion at 10 kpc. The case of inverted hierarchy 
and large 6*13 is considered. The results are obtained for 10 
(solid), 15 (dashed), 19 (dash-dotted) and 29 (dot-dot-dashed 
line) MeV positron energies. 



Chooz limit. They correspond to the case of the adia- 
batic conversion at the high density resonance [24| . The 
other value is chosen as an example of the non-adiabatic 
regime. Note that sin^26'i3 = 10^^ — 10^^ arc the small- 
est values that can be reached in accelerator experiments 
even with very long-term projects (25| . 

For the neutrino luminosities at the neutrino-sphere 
we make the usual assumption of equipartition of energies 
among all neutrino flavors and that they decrease as = 



L^g X exp(— i/r) with = 10 erg • s~ and r = 
3.5 s. Equal luminosities are appropriate for the cooling 
phase of the neutrino signal upon which we are focusing; 
during the accretion phase the and De luminosities are 
substantially brighter than the [l^. We consider the 
hierarchy i.e. (i?^^) < (Ejy^) < {E^^) with typical values 
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of 12, 15 and 18 MeV respectively. 

Figure [1] shows the evolution of probabilities for an 
anti-neutrino that is produced at the neutrino-sphere as 
an electron- type. One recognizes the synchronization re- 
gion and the rapid bi-polar oscillations near the neutrino- 
sphere. The electron anti-neutrino probability becomes 
very small implying that due to the vv interaction, the 
corresponding fluxes have swapped with the v^^Vt neu- 
trino fluxes and has turned "hot" at this point in their 
propagation. 

To help in understanding the results, let us now con- 
sider a supernova explosion located at 10 kpc from Earth. 
At the early times (t < 1 s) , for the inverted hierarchy and 
large ^13, anti- neutrinos undergo an adiabatic MSW res- 
onance and have a "cold" spectrum on Earth^ (Figure[2]). 
When the shock wave passes through the MSW high res- 
onance region, important modifications of the neutrino 
fluxes occur. The presence of the shock wave renders 
the neutrino flavor conversion non-adiabatic. Therefore 
the neutrino spectra on Earth remains "hot" (Figure [2]). 
Note that the adiabatic and non-adiabatic spectra will 
cross at some energy, which in this case, is E^, ~ 20 MeV. 

Our predictions for the positron time signal associated 
to inverse beta-decay in a detector are shown in Figure 
[3l Let us consider the case of inverted hierarchy and 
large third neutrino mixing angle. At early times, the 
conversion is adiabatic and neutrinos with less than 20 
MeV will produce a number of positrons determined by 
the "cold" spectrum. This number will decrease when 
the shock wave renders the resonance conversion non- 
adiabatic and the neutrinos reach the Earth with a "hot" 
spectrum. This will show up in the positron time signal 
as a dip. Its depth is energy dependent. (Note that if 
the w interaction is absent, the time signal would show 
a bump instead, as discussed in [l^.) For energies larger 
than 20 MeV, since the relative number flux of "hot" and 
"cold" spectra interchange, the dip turns into a bump. 

Figure[4]presents our prediction for Super-Kamiokande 
(22.5 ktons) for which about ten thousand events are ex- 
pected. If one takes 0.5 s time bins, and a lower energy 
bin between 10 and 19 MeV and a higher one, above 25 
MeV, one can see that a bump (dip) can be distinguished 
by the pure exponential behaviour at 3.5 (1) sigma for 
the upper (lower) energy range. Obviously such a signa- 
ture could be observed with exquisit precision if a next- 
generation large-size detector such as MEMPHYS (UNO, 
Hypcr-K) is considered [S^] . The large number of events 
would allow to use small energy bins and follow smoothly 
the transition from bump to dip. 

Finally, let us discuss the sensitivity to the hierarchy 
and upon the third neutrino mixing angle (Figure [5]). If 
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FIG. 4: Positron time signal in Super-Kamiokande, taking 0.5 
s time bins, and two energy bins : 10 to 19 MeV (upper fig- 
ure) and above 25 MeV (lower figure) . An exponential decay 
calibrated on the first two time bins is shown for comparison. 



the hierarchy is inverted and ^13 is large the positron 
signal as a function of time presents a dip as discussed 
above. If O13 is small, no flavor conversion occurs in 
the high density resonance region and the electron anti- 
neutrino fluxes remain "hot" . This leads to simple expo- 
nential decay, in the positron time signal. On the other 
hand if the hierarchy is normal the anti-neutrino fluxes do 
not change traversing the whole star and remain "cold" 
during the explosion. In this case the positron time sig- 
nal again is an exponential decay but more positrons are 
produced compared to the inverted hierarchy case. 
Conclusions We have performed the first calculation to 
include the most recent developments of neutrino propa- 
gation in dense media such as core-collapse supernovae. 
Our numerical results include the neutrino-neutrino in- 
teraction on one hand and use evolving density profiles 
from realistic simulations which include multiple reso- 
nances on the other. Our results show that the inter- 
play between the neutrino-neutrino interaction and the 
shock wave can, indeed, be understood in terms of spec- 
tral swapping. We have explored the positron time signal 
related to electron anti-neutrino on proton scattering, the 
dominant detection channel in Cherenkov and scintilla- 
tor detectors. We have shown that the positron event 
rate possesses a characteristic time signal that depends 
upon the neutrino hierarchy and third neutrino mixing 
angle. For an inverted hierarchy and a large third neu- 
trino mixing angle the event rate is found to decrease (in- 
crease) midway through the supernova neutrino signal for 
low (high) neutrino energies. Such a signature is observ- 
able in current and future large-size neutrino detectors 
at present under study. Note that this general predic- 
tion agrees qualitatively with the observed gap of low 
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^ Note that the spectra mix shghtly due to the 612 rotation at the 
L resonance. 



energy events in the supernova 1987A data 
although, with so few observations, an emission model 
with no shock effects is marginally compatible [3l| . The 
signature relies upon: a) that the proto-neutron star is 
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FIG. 5: Positron time signal associated to inverse beta-decay 
in a detector, for a 12 MeV positron. The results correspond 
to a normal hierarchy (solid) or an inverted hierarchy with 
large (dot-dashed) or small (dashed) third neutrino mixing 
angle. The presence of a dip is typical of the whole range 
sin^ 6^13 = 0.1 - 10"". 



brighter in Df. and Dx at low energies, b) that the hier- 
archy is inverted so that collective effects swap the j>e 
and Vx spectra prior to the H resonance, c) that neutrino 
propagation through the progenitor profile is adiabatic, 
d) the shock reaches the H resonances while the super- 
nova is still luminous. If, for example, the luminosities 
of De than Vx at low energy were comparable then there 
would be little decrease in the positron event rate mid- 
way through the signal but, on the flip side, a greater dif- 
ference in luminosity exagerates the event rate decrease. 
This signature might be robust in the presence of tur- 
bulence, although further investigation is required. The 
reasoning is that the density profiles used here can be 
thought of as equivalent to the 'average' profiles of Fogli 
et al. [s^. When these authors added turbulence to the 
post shock region of their profiles they found that the 
size of the shock effects were muted but, more impor- 
tantly, they were not removed. Finally the observation 
of the cross-over energy provides valuable information 
above the original neutrino fluxes. By measuring differ- 
ent positron energies in a possible future galactic core- 
collapse supernova explosion one might learn if the third 
neutrino mixing angle is within (smaller) the window of 
achievability of present (future) terrestrial experiments. 
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